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The units tell you what to do 

Simon Brown 


Abstract: Many students have some difficulty with calculations. Simple dimensional analysis 
provides a systematic means of checking for errors and inconsistencies and for developing both 
new insight and new relationships between variables. Teaching dimensional analysis at even the 
most basic level strengthens the insight and confidence of students, and provides them with a tool 
that can be applied in almost any context. 
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1. Introduction 

Every science teacher has been confronted with a student unable to work out how to do certain 
calculations. Sometimes it is a simple problem, like calculating the number of moles in a given 
volume of solution of a specified concentration, but it can be more complicated. The decline in the 
mathematical skills [1, 33] and insight of students has increased the frequency with which this 
problem is encountered. Given this, strategies for addressing such problems have become increasingly 
important. One approach is to use the units as a guide. 

I often repeat the title of this paper to my students in the hope that they will remember this idea if 
nothing else. The assertion is based on the concepts of dimensional analysis (DA), a technique often 
taught to chemists, physicists and engineers, but less frequently to biologists. At its simplest, DA can 
be used to check for errors in calculations or equations, but with very little more sophistication it can 
also enable one to work out the general form of an equation rather than having to remember it or look 
it up. The most sophisticated application facilitates the fonnulation of new relationships between 
variables: for example the relationship between particle size and the wavelength of light in Rayleigh 
scattering was determined using DA [32]; and it has been used by many other people, including 
Maxwell [24], Fourier [16] and Einstein [14], 

If the correct use of units can be infonnative, an incorrect use can be catastrophic. Perhaps the best 
known example of the latter is the Mars Climate Orb iter Tost’ by NASA in 1999 because two groups 
of engineers used different units preventing correct data transfer [21], Aviation is especially rich in 
such confusion, for example in 1983 an Air Canda flight ran out of fuel on the way to Edmonton 
because a non-metric fuel density (1.77 lb L' 1 ) was used to convert a mass of fuel in kilogrammes to a 
volume in litres [2]; in 1997 a Korean Air flight flew into a hillside as it approached the airport at 
Guam, apparently because the flight crew were unclear whether they were working in feet or metres 
[27]. Medical professionals sometimes have difficulty calculating drug doses because they have to 
interconvert units leading to the accidental administration of incorrect doses of medication [36, 37]. 
On at least two occasions American long and triple jumpers have encountered problems in 
international competition because metric units are used, rather than the imperial units (feet and inches) 
used in the United States. For example, Melvin Lister did not reach the triple jump final at the 2004 
Olympics. Later he said “They made me use a metric tape for my run-up and I'm used to imperial, ft 
ruined my trust.” [29] and, perhaps more understandably, “Nobody told me they were only going to 
have metric out there. I couldn't figure out what my mark was.” [12], 

Here I outline the basis of DA and provide some examples of its application. In the first instance, I 
give some very simple illustrations of the fundamental ideas needed to develop and analyse technical 


Received 31 October 2009, accepted 21 December 2009, published 31 December 2009 






92 


Simon Brown 


problems (sections 2 and 3). Finally, I outline a more fonnal approach and give some more examples 
(section 4), and then describe the value added from teaching DA (section 5). 

2. The fundamentals of dimensional analysis 

The general rule is very simple: the dimensions of all parameters related by addition, subtraction and 
an equals sign have to agree. Specifically, the dimensions on either side of ‘=’ must be the same 
(obviously, it makes no sense to write something like 1 mL = 1 mm); if numbers are added or 
subtracted they must have the same dimensions (for example, it makes no sense to write 5 mL + 
2 mm); and if numbers are multiplied (or divided) the dimensions are just combined according to the 
mathematical operations (2 mm x 2 mm = 4 mm 2 or 20 km h' 1 / 2 km = 1 h' 1 ). Moreover, it is not 
possible to apply some mathematical functions to a unit (for example, it is not reasonable to write 
log(l mm), sin(l mm) or exp(l mm)) so the arguments of those functions have to be dimensionless. A 
useful practice is to replace all SI prefixes with the appropriate multiplier (for example 1.5 mm = 1.5 x 
10' 3 m or 1 ML =1 x 10 6 L), which not only maintains dimensional consistency, but provides an 
explicit connection with ‘scientific notation’. 

Maxwell [24] introduced the convention that the dimensions of a variable ( x ) are indicated using 
square brackets ([*]). In chemistry such square brackets usually denote concentration, but I will use 
Maxwell’s convention here and write the concentration of x as c x . 

Naturally, dimensional consistency and dimension interconversion can be analysed using this idea 
[34], but one can also use it to work out how to combine variables correctly. For example, the 
common unit of concentration (c) is molarity ([c] = mol L' 1 ), but molality (in) is also used ([m\ = mol 
kg" 1 ) and these units are connected by the density of the solvent ([/)] = kg L" 1 ). To convert a molarity 
to a molality assume that 

[m] = [c] a x lof 


or 

mol kg" 1 = (mol L" 1 )" x (kg L' 1 )* = mol 17 * L “ * * kg' 7 
and, by equating the exponents of the units on each side 

mol 1 = mol" => a = 1 and kg" 1 = kg* => b = —1 

but, obviously, since L does not appear on the left L a b = L°, so -a b = 0, consistent with these 
values of a and b. Of course, this means that m = c x p . While this is a trivial example, the same 
approach works with much more complicated problems and I shall make use of it below. 

The same approach works equally well if the dimensions are reduced to their essential components: 
length (L), mass (M), time (T) and temperature (0). More complex units are just combinations of 
these fundamental dimensions, for example 

1. speed is just distance (L) moved in a given time (T), so [speed] = L T' 1 ; 

2. density is just a mass (M) in a given volume (L 3 ), so [density] = M L" 3 ; 

3. energy is measured in joules (J) and 1 J = 1 kg nr s' 2 (the energy needed to accelerate 1 kg at 
1 m s" 2 over 1 m), so [energy] = M L 2 T" 2 ; and 

4. work (measured in watts, W) is the rate at which energy is expended, so [work] = [energy] T" 1 = 
M I. 2 '!" 3 . 

Not having to worry about all the details of the specific units can make it easier to work out what 
calculation is actually required. 

3. Devising new relationships and gaining insight 

A simple application of DA can lead to apparently surprising insights that challenge students to 
question their assumptions, leading to deeper understanding. For example, one often sees 
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pH = -log 10 (c H+ ) 


( 1 ) 


or, equivalently, 


c H+ =10' pH , (2) 

relating pH to the concentration of protons (c H +). By convention [c H +] = mole L' 1 , rather than any of 
the other concentration units that might be employed, and, since [pH] = 1 and therefore [10" pH ] = 1, (2) 
is clearly dimensionally inconsistent (as is (1), since (1) and (2) are equivalent). As the argument of 
the logarithm must be dimensionless (because log(mole L" 1 ) is a nonsense, as Sorenson [31] 
recognised), there must be a factor that converts c H + into a dimensionless equivalent (« H +, say), so that 
(1) is actually 


pH = -logio(a H +) • (3) 

In practice, a H + is the thermodynamic activity of protons, to which c H + is related by the apparent 
activity coefficient (Jh+'), so u H + = Yn+'cn+ and [yn+] = L mol" 1 . It is more usual [7] to assume a 
dimensionless activity coefficient y H + = cqYh+', where Co is a standard concentration (= 1 mol L" 1 ), so (1) 
becomes 


pH —log 10 (a H +) —logio(^+'c H +) —logio(XH+CH+/co) —logio(TH+) - logio(cH+/co). 

Only in an ideal solution does yk+ = 1, in which case 

pH = —logio(c H +/c 0 ), 

so values of ;/n- ^ I indicate deviation from ideality. The realisation that (1) is dimensionally 
inconsistent, and the resolution of that inconsistency, prompts students to question and helps them to 
develop considerable insight. The simple reasoning leads to a deepening of student understanding. 

In addition to dimensional inconsistency, there are at least three other problems with (1) and (2). First, 
there need be no unique relationship between c H + and pH, because Yh+ need not be constant (for 
example the mean activity coefficients of various acids depend on the concentration of the acid [19]). 
McCarty and Vitz [25] give several examples for which (1) is incorrect in this respect, and suggest that 
students should be taught (3) rather than (1), but concentration is a challenge to many students [17] 
and activity is even more difficult. The second question is the significance of the ‘p’ in pH. If pH is 
related to the concentration of H + , then what has the ‘p’ got to do with it and why does it mean that pH 
decreases as c H + increases? It is of little help to most students to relate the ‘p’ to potenz, puissance or 
power [22, 28] and suggesting that ‘p’ is related to ‘-logio’ is often only marginally better and reduces 
the whole issue to one of calculation. The third problem, related to the second, is that the significance 
of ‘logio’ is often unclear to students because of their relatively poor mathematical literacy. In many 
classrooms the use of pH has been reduced to the question of which button to press on a calculator 
[35]. Whether (1) or (3) is taught to students, this fundamental difficulty remains and the question of 
whether the argument of the logarithm should be c H + or a H + is actually secondary to the more 
important consideration of mathematical competence. 

A simple, and well known, extension of this example relates to the dissociation of an acid (HA) 

HA - H + + A”, (4) 

for which the mass action ratio (Q) is usually written 

Q = Ch+CaVcha, (5) 

where ca- and c H a are the concentrations of A~ and HA, respectively, and the equilibrium constant is 
K a . At equilibrium, Q = K a and one often sees this written as the Henderson-Hasselbalch equation [11, 
30] 

pH = p K a + logio(c A -/cHA), (6) 

where p K a = log| 0 (A„) and pH is given by (1). Since [pH] = 1 and because the argument of the 
logarithm is also dimensionless, [p/T,] = 1. Consistent with the argument leading to (3), [A"„] = 1 and 
so [ Q\ = 1 [10], but it appears from (5) that [Q\ = mol L" 1 , indeed K a values are often reported as 


Volume 2, Number 4, 2009 



94 


Simon Brown 


having some specific dimension [9], despite recommendations to the contrary [10]. For an electrolyte 
that dissociates into v+ cations and v anions, the activity is defined as a so i u t e = a ' a v = a± v , where v 
= v+ + v . Since the activities of individual ions can not be measured, the (geometric) mean activity 
coefficient y± = (y+y_) Vv is used. For the acid dissociation reaction (4), v+ = = 1 and y+ = (y a+ y A ) v \ 

which varies with the concentration of the acid (Figure 1). Based on this, the Henderson-FIasselbalch 
equation (6) becomes 

pH = p K a + logio(c A VcHA) + logio(r±) 5 

where the third term on the right can be thought of as representing deviation from the behaviour of an 
ideal solution. 



m (mol kg' 1 ) 


Figure 1. Activity coefficients (y±) of some hydrogen halides (HCl, HI, HBr, HF) [19] as a function of molality 
(m) and the limiting (7) and extended (8) Debye-Hiickel approximations (dashed and solid curves, respectively). 


For molalities ( m ) of less than about 0.05 mol kg" 1 (Figure 1), y± can be estimated using the Debye- 
Hiickel limiting equation 

lny ± = -|z + z_|«V7, (7) 

where z+ and z_ are the valences of the cation and anion, respectively, (so z+ = 1 and z_ = -1 for the 
dissociation of an acid) and 


a = yjlxN A ps 3 , 


where N A is Avogadro’s constant, p is the density of the solvent and 

2 


47r£ 0 £ r kT 


is a distance characteristic of the solvent which depends on the dielectric constant ( e r = 78.54 for 
water), the charge of an electron (e = 1.6021 x 10" 19 C), Boltzmann’s constant (k = 1.38 x 10' 23 J K" 1 ), 
the absolute temperature (T= 298.15 K) and the permittivity of free space (i' 0 = 8.85 x 10" 12 C V" 1 m" 1 ). 
Using these values, s = 0.7136 mn and a = 1.171 (kg mol" 1 ) /i . At higher molalities (m < 0.1 mol kg" 1 , 
the extended Debye-Hiickel equation [13] 


In r± = 


|z + z_| cx-JJ 
1 + 2 a(a/s)s/J 


( 8 ) 


Acta Didactica Napocensia, ISSN 2065-1430 










The units tell you what to do 


95 


is more reliable (Figure 1). In (8), a is the average smallest distance between two ions (a = 0.30-0.35 
mn [18]). Both (7) and (8) depend on the ionic strength (/) of the solution 

/ =0.5 ^rrijZj 

j 

[23] where ntj is the molality of the ion j which has valence z, and so [7] = mol kg' 1 . It is clear from 
Figure 1 that neither (7) nor (8) is adequate at high concentrations, and also that they do not work for 
hydrofluoric acid (HF), which behaves quite differently from other hydrogen halides (HC1, HBr and 
HI) because of the solvation of F~ by HF at even quite low concentrations [26]. Indeed, whereas HC1, 
HBr and HI are strong acids (p K a = -7, -9 and -11, respectively), HF is a weak acid (p K a = 3.45) at 
low concentrations. 

The usual means of deriving an expression for the activity coefficient is based on the observation that 
for a solute i the chemical potential (//,) can be viewed as the sum of the chemical potential of the ideal 
solute (//,-,ideal) and a term representing non-ideality (//,-,non-idcai) which is 

/7.non-ideal — RT Iny,-, 


where R is the gas constant. If the non-ideality arises from electrostatic interactions, then 


In Yi=~ 


2 2 r, 

ZjZ P 

%ne^e r kT 


( 9 ) 


where [:> is the Debye length, which is a measure of the approximate extent of the electrostatic field 
of the ion, 


nl e ST 2 
13 


H z i 

and to convert [«,] = m' 3 to [/«,] = mol kg' 1 we use m, = N A pcJ{ 1CT 3 nT L' 1 ) to obtain 


„2 A »pe 2 2iVape 

P~ = . ^- 1 ■ 


10 3 £ 0 £,.kT‘ 


10 3 £ 0 £ r kT 


Substituting this into (9) yields 


In Yi = 


2 ttN a P jj 

{4K£ 0 £,.kT)P 2 V 10“ 3 


= —zf ayfJ , 


but, for the acid dissociation reaction (4), = (j+y )'' 2 or In y+ = '/i(ln y+ + In ;/) and z+ = 1 and z = I. 

so the Debye-Huckel limiting equation (7) is 

In y ± =-a4l (10) 

which, in water at 298 K, is In = 1.174 f 2 or logi 0 y+ = 0.509 1 /2 (Figure 1). 

Even this derivation is relatively demanding, but an approximation of the Debye-Huckel limiting 
equation can be obtained quite simply using DA, as will be demonstrated in the next section. 


4. A more formal version of dimensional analysis 

Buckingham [8] described a more fonnal approach to DA that can be stated in a few steps: (1) count 
the number of variables (N v ) and the number of independent dimensions (N d ); (2) fonn N r - N d 
dimensionless groups ( 0„ i = 1, 2, ..., N v - A)/); and (3) the most general relationship between the 
variables is 

0 1 =F(0 i ,i = 2, 3, ...,N V -N d ), 

where F is a dimensionless function. Each dimensionless group is fonned by multiplying the variables 
involved (say x and y) raised to separate powers (say a and b, as illustrated in section 2) to form 0 = 
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x a y h such that [ 0\= I. There may be several possible ways of forming the ft, and one can choose from 
the options. 

In the case of the Debye-Hiickel limiting expression (7), if In y± depends on electrostatic interactions 
between ions, then it involves e, s 0 , s,, k, T, p, N A , and 7, for which 

[In y ± \ = 1, [e] = Q, [£<,£,] = Q 2 (M L 3 T 2 ) 4 , [k] = M L 2 T 2 © 4 , 

[7] = 0, [p] = M L- 3 , [N a ] = n 1 , [7] = n M" 1 

and N v = 8 and N d = 6, so two dimensionless group can be formed. The first is just ft = In y+ and the 
other is 


[ft] = 1 = Q a x (Q 2 (M L 3 T" 2 ) 4 ) 6 x (M L 2 T' 2 0 4 ) c x0^(M L' 3 / x ( n 4 ) g x ( n M' 1 )*. 
Collecting tenns yields 

Q fl + 2i X M - h + c+f - h X L" 36 + 2c - 3/ x j2b-2c x Q-c+d x n -g+h = J 

but since all of the exponents must equal zero, then 


"1 2 0 0 0 

0-1 10 1 

0 -3 2 0 -3 

0 2 -2 0 0 

0 0-11 0 

v 0 0 0 0 0 


0 

0 

0 

0 

0 

-1 


-1 

0 

0 

0 

1 


A 

V 


b 


c 


d 


f 


g 

/ 

J 1 a 


= 0 . 


This system of equations can not be solved explicitly for all the variables because there are fewer 
expressions than variables, but it is clear that a = -2b, b = c = d and/= g = h. Taking/= g = h = 1, 
then b = c = d = - 3 and a = 6 from which 


0 7 = 


r 2 \ 

e 

e 0 e,.kT 


N aP I 


and the general form of the expression is 

\ny ± = F 


r f 2 V ^ 

e 


y £ 0 £ 1 -kT y 


N aP I 


( 11 ) 


which is clearly similar to (10), but it required no complex mathematics. If it had been assumed that/ 
= g = h = 0.5, then the Debye-Hiickel limiting equation (7) would have been approximated even more 
closely. 

As a final, slightly more complex illustration of the power of DA, consider Bjerrum’s [4] expression 
for K a . Bjerrum [4] combined the Debye-Hiickel equation (7) with the notion of mass action (6) and 
obtained an expression for the equilibrium constant based on the electrostatic interaction between ions 
a distance r apart and a distance (7?) representing association of ions. Assuming that the variables 
involved are 

[R] = L, [r] = L, [K a ] = 1, [e] = Q, [£b£>] = Q 2 (M L 3 T 2 ) 4 , 

[k] = M L 2 T" 2 0 4 , [7] = 0, [p\ = M L" 3 , [/V,] = n 4 , [7] = n M 4 
then N v = 10 and N d = 6, so four dimensionless groups can be formed 

ft = Ka, ft = N A Ipr 3 , ft = e 2 /r£ 0 £,kT, ft = r/R. 

Clearly ft and ft are linked, because (11) can be written In y± = F(ft ft 3 ), so ft = ft 3 F'(ln y ± ) and 
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K a = G(N A Ipr\ e 2 /r£ 0 £,kT, r/R ) = G(F '(In r± )(r£ 0 £,kT/e 2 )\ e 2 /r£ 0 £ r kT, r/R). 

Using these dimensionless groups, Bjerrum’s [4] expression can be written as 

, 4^ , 

K =-ft ft 3 f y 4 exp (v)dy , 

1000 / J a 

so even without solving the integral it is clear that K a = G{ ft, ft, ft), as shown using DA. 

These two examples illustrate the extraordinary power of DA to provide insight while using only 
simple algebra. Most significantly for a teacher, it facilitates the analysis of a problem from a 
scientific perspective without requiring the use of complex mathematics [5, 20], at least in the first 
instance. While DA will not always yield an answer, it can be relied on to generate insight into the 
problem. This is the basis of the assertion that the units tell you what to do. 

5. The value of teaching dimensional analysis 

Teaching even the simplest level of DA is powerful, because it develops (i) questioning; (ii) the habit 
of checking one’s own work, particularly in the sense of trying to assess whether an answer is 
reasonable; (iii) skills in the use of the SI prefixes and the link with scientific notation, neither of 
which is well understood by many students; (iv) fundamental algebraic skills; and (v) a powerful 
technique that can be applied in any discipline. A simple application of DA to familiar concepts, such 
as pH (section 3), elicits useful discussions and motivates students to think differently about what they 
think they know. Ultimately, the habits acquired by using DA contribute to the metacognitive 
development of students. Even if students do not appreciate DA initially, they are eventually 
persuaded of its value. 

6. Conclusion 

Dimensional analysis is a valuable tool that should be taught much more widely. Even used in a very 
simple way it is powerful and useful. Many examples of the use of DA can be found in the literature 
[3, 5, 6, 15, 20], but with a little practice it is not difficult to develop new examples. As with much in 
science, the ‘secret’ is allow yourself to think simply and let the units tell you what to do. 
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